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We propose an efficient mechanism for the evaporative cooling of trapped fermions directly into 
quantum degeneracy. Our idea is based on an electric field induced elastic interaction between 
trapped atoms in spin symmetric states. We discuss some novel general features of fermionic evap- 
orative cooling and present numerical studies demonstrating the feasibility for the cooling of alkali 
metal fermionic species 6 Li, 40 K, and 82,84 ' 86 Rb. We also discuss the sympathetic cooling of fermionic 
O ' hyperfine spin mixtures, including the effects of anisotropic interactions. 
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I. INTRODUCTION 



The success of atomic Bose-Einstein condensation 
(BEC) has opened a major research area |J-§] in weakly 
interacting quantum gases B. Experimental break- 
throughs have been achieved thanks to several remark- 
able advances in atomic, molecular, and optical physics. 
In particular, the laser cooling and trapping of neutral 
particles |^](| , and the ingenious application of the evap- 
orative cooling technique 
^i . Evaporative cooling was first developed for magneti- 
1 cally confined atomic H and is based on a simple 
principle seen in everyday life In a quasi equilib- 

C***'* ■ rium system the particles most likely to evaporate are the 
most energetic. In doing so, the particle leaves behind a 
system with a lower average energy per particle, i.e., a 
cooled system. Less appreciated perhaps, is that forced 
evaporative cooling is the key to reaching quantum de- 

Sgeneracy in alkali metal systems. A sufficient condition 
for successful evaporative cooling is the threshold "run- 
away" requirement, when an increased collision rate is 
achieved despite the loss of atoms. This technique was 
adapted to alkali atomic systems in the first generation 
of BEC experiments |l(| . It has since been demonstrated 
to work over a wide range of temperatures and densities, 
and has been applied to all successful BEC experiments 
this far (J. 

With such a powerful technique one may ask why it 
has not yet been used to create degenerate atomic Fermi 
gases ? A careful examination of current magnetic trap- 
ping set-ups helps to clarify the fundamental problem. 
Alkali metal atoms have two degenerate manifolds of hy- 
perfine ground states / = J± 1/2, due to the coupling 
between the nuclear / and electronic spin s (s — 1/2). A 
typical magnetic trap arrangement confines the low field 
seeking Zeeman state with maximum spin polarization 
rrif — f (at the moment we assume a single rrif value 
is trapped, mixtures of two different mf values will be 
discussed later). Eigenstates of F 2 , where F = f + f is 
the total hyperfine spin of two colliding atoms, have def- 
inite spin exchange symmetry given by (— l) 2 f~ F = +1. 



Thus only the partial waves for which the relative or- 
bital angular momentum I is odd contribute to the an- 
tisymmetrized scattering amplitude. From low energy 
atomic collision theory, we know that the partial wave 
phase shifts are 6i(k) oc k 2l+1 for a short ranged poten- 
tial, typical of alkali metal atoms in their ground state. 
Therefore only the s-wave (I = 0) has a non-zero con- 
tribution in the low energy (k — > 0) limit, characterized 
by a scattering length a sc = — lim/^o Sa(k)/k. One then 
obtains the elastic cross section ob = 87rOg C for bosons, 
while for fermions op — * 0, in the limit k — > 0. Without 
the thermal equilibration produced by elastic collisions, 
evaporative cooling can not be directly applied to a single 
fermionic species. This is the same fundamental reason 
that BCS states for trapped fermionic gases are not easily 
accessible ■ Obviously, if two or more spin compo- 

nents are simultaneously trapped, spin-statistics does not 
prohibit s-wave collisions between the spins at low energy, 
so that sympathetic cooling is in principle possible. 

In this paper we discuss an efficient mechanism for the 
evaporative cooling of trapped fermions. The mechanism 
employs an external electric field to polarize atoms, en- 
abling cooling at low energy. We also discuss the sympa- 
thetic cooling of fermionic mixtures. The paper is orga- 
nized as follows. In section II we begin with a discussion 
of the major results for collisions between atoms in an 
external electric field |E| . We then outline our proposal 
for the evaporation kinetics of fermions, closely following 
the approach developed in for bosons. We present 

our numerical studies and estimate the time scales for 
fermionic species 6 Li and 40 K. In section III we discuss 
sympathetic cooling of fermionic spin mixtures, along the 
lines of the Jila experiment |L6) . Finally, in section IV 
we present our conclusions. 



II. COLLISION KINETICS OF ATOMS IN AN 
ELECTRIC FIELD 

As shown in detail in ref. (l3|| , polarized atoms in an 
external electric field interact through the dipole interac- 
tion 

Cw 

V E (R) = -^P 2 (co S 6), (I) 

in addition to the usual interatomic potentials. Here, 
Ce = 2E 2 a^ (0)af (0) is the electric field induced dipole 
interaction coefficient and af' B (0) are the static atomic 
dipole polarizabilities of the two atoms denoted by A and 
B, respectively. The angle between the directions of the 
electric field ( of amplitude E ) and the internuclear axis 
is denoted by 6, and P2 is the Legendre polynomial of or- 
der 2. The presence of this anisotropic interaction term 
completely changes the low energy collision physics. In 
particular all partial waves now contribute to the low 
energy scattering cross sections since one can prove ana- 
lytically that 5i(k) ~ k for all I Jf3|. The unsymmetrized 
low energy scattering amplitude is 



(2) 



Iml' 



where t\™ are the reduced T-matrix elements, the multi- 
channel analogues of the scattering length. Here the in- 
cident relative momentum of the two colliding atoms is 
denoted by k = kk, and k' = kR is the relative momen- 
tum after the collision. 

For a gas of atoms in an electic field we can now derive 
the collision integral for the quantum Boltzman equation 
(QBE). The QBE is given by 



d P „ 

at m 



V ? C/(f) • Vp f(f, p, t) = l(r, p, t), (3) 



where f(r,p,t) is the phase space distribution function, 
and U (r) the external trapping potential. Using Eq. (Q) 
we arrive at the following collision integral 
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[l + Vf(f,PiW + vf(r,P2)]f(r,p 3 )f(r,M) , (4) 



where is the symmetrized, unsymmetrized, or the 
anti-symmetrized scattering amplitude for bosonic (77 = 
1), classical (77 = 0), or fermionic atoms (77 = —1), re- 
spectively. We have used the following notation. Two 
atoms with incoming momentum p 3 and p 4 collide and 
produce an outgoing pair with momentum pi and p 2 . The 
scattering amplitude from an initial relative momentum 
q = (p 3 — pi)/ 2 to a final momentum q' = (pi — p%)/2 
is given by F(q,q'). For elastic collisions one has q = q' 
in addition to the conservation of the total momentum 
P = P3+PA = P1+P2 [assuming the collision to be local- 
ized in space at a fixed U(r)] and energy 63 + 64 = 61+62. 
The differential cross section for scattering into the (rel- 
ative) solid angle dq' is given by da = \F(q,q')\ 2 dq', 
and \F(q, q')\ 2 = \F(q, q')\ 2 for the low energy collisions 
under consideration. The approximation of localized col- 
lisions within a small spatial region compared with the 
variations of U(r) allows us to neglect the P dependence 
of F. 

Assuming an ergodic distribution for the system, we 
introduce the energy dependent single particle distribu- 
tion function n(e, t) and obtain the kinetic equation in 
ergodic approximation p"4| , 

p(e 4 )-^-n(e 4 , t) sa / de 4 de 2 de 3 5(ei + e 2 - e 3 - e 4 )p(min[ei, e 2 , e 3 , e 4 ]) 

at 7r 2 7i J 

n(ei)n(e a )[l + vn{e 3 )][l + r]n{e 4 )] - [1 + »?n(ei)][l + f]n{e 2 ))n{e 3 )n{e 4 ) , (5) 

where p(e) is the density of states. Our result is for- 
mally similar to the standard result for isotropic s-wave 
collisions in the absence of an electric field |I^,|l5|, when 
the s-wave collision cross section is replaced by the effec- 
tive scattering cross sections o^. In reaching the form of 
Eq. (El), the collision integral has been approximated 



by assuming that the fast varying part in the cos 9^ 
and cos#g-/ integration comes from the delta functions 
6[e — pf /2m — U (r)] rather than the scattering amplitude 
F ( j )1 thus we have extracted the scattering cross section 
factor out of the integrand, and replaced it by its average 

1,1 even/odd 

J dq'dq-\F{qA')\ 2 = 2 x 4tt £ |t^'| 2 = *„, (6) 

where the additional factor of 2 occurs in the Bose/Fermi 
case due to the symmetrization/anti-symmetrization. A 
specific discussion of p-wave scattering is presented later. 

For highly confined atoms inside magnetic traps, cor- 
rections due to discrete level structure need to be in- 
cluded. When a long range type of potential ([[]) is in- 
volved, one can not in general approximate the atom- 
atom interaction potential by a contact 5{f — ?) form as 
is usually adopted in the s-wave collision models. Spe- 
cial care is needed to find the scattering matrix clement 
for atoms in different trap states. The quantum kinetic 
equation now takes the form (for bosons and fermions) 

dn i x 

l,k,l 

(njn^l ± n ? )(l ± nj) - njn y (l ± n r )(l ± n g )) , (7) 

where — n(e^) is the population in level = h(i x Lu x + 
iyijJy + i z uj z ) with ui x , uj y , and lo z the trap frequencies in 
the x, y, and z directions. 

For s-wave collisions, when the contact potential is 
valid, the collisional matrix element is given by 



7-? t r oc 

! i,j;k,l 



$* (r) $* (r ) $ t (f) $ Ar) dr 



(8) 



where denotes the wave function for trap state i = 

(i x ,iy,iz)- For the anisotropic collisions discussed here 
we have 



Im.l'm 



(9) 



_ 2tt 1 
with 

4j'(P)= j dq'^(p/2 + q r')^/2-?)Yt m ,{q% (10) 

where is the wavefunction in momentum space. 

Clearly, apart from final state quantum statistical fac- 
tors, the collision kernel on the rhs of Eqs. (||) and (0) 
takes the same form as that for classical atoms with s- 
wave collisions. For fermions, the effective cross section 
<7_ provides an opportunity for investigation of the er- 
god ic kinetics as has been done previously for bosons 
Therefore the net effect of the electric field is 
to polarize atoms, modify their collisional properties and 



thereby induce a non-zero elastic cross section at low en- 
ergy. The electric field dependence of the effective cross 
sections a = <r± allows one to control the time scale of 
evaporation by adjusting the strength of the electric field 
E. In addition, since the microscopic collision processes 
for polarized atoms are anisotropic, kinetic motions in 
different spatial directions are mixed. The validity of the 
ergodic approximation is thus reinforced. 

A. Results 

In this subsection we report our numerical studies of 
fcrmionic evaporative cooling in an electric field. We 
focus primarily on the qualitative differences between 
fermions and bosons in this paper, and the fundamental 
issues relating to fermionic evaporation. We will approx- 
imate the scattering matrix element 7jj.g f by the equiv- 
alent s-wave contribution, retaining only the leading £qq 
contribution in (^). We believe that this approximation 
should not produce any qualitative differences from the 
full model in the results presented. 

In the early stages of evaporative cooling, when the ef- 
fect of phase space blocking term [1 — f(r,p)] or (1 — n e ) 
is unimportant, the kinetic term reduces precisely to that 
of a classical gas with an effective cross section rr_ . The 
evaporative kinetics of such a case have been extensively 
studied . The classical run-away condi- 

tions can be achieved provided loss mechanisms are con- 
trolled. For 6 Li, 40 K, and 82 ' 84,86 Rb, we have performed 
detailed multichannel collision calculations for atoms in 
an electric field. For fields of the order of a few (MV/cm) 
we find that one can readily induce an effective scatter- 
ing cross section as large as for their bosonic counter- 
parts 7 Li, 39 K, 41 K, 85 Rb, and 87 Rb. Our results are 
displayed in Figs. |l|, |[ and ||[ The case of 6 Li is particu- 
larly interesting, since the triplet scattering cross section 
reaches the field-free value of bosonic 7 Li (87ra|,) with 
a T = -27.6 (a.u.) at a value of E ~ 3 (MV/cm) [§. 
For 40 K, it reaches the field free bosonic isotope value of 
a T ~ 80 (a.u.) for 39 K and a T ~ 300 (a.u.) gl for 41 K 
at values of E ~ 2 (MV/cm) and E ~ 1 (MV/cm), re- 
spectively For 82 ' 84,86 Rb, it reaches the field free bosonic 
isotope value of or ~ 369 (a.u.) for 85 Rb and ar ~ 106 
(a.u.) for 87 Rb @ at values of E ~ 0.8 (MV/cm) and 
E ~ 1.5 (MV/cm), respectively. 

The promising scenario discussed above can not be 
simply extrapolated to the regime of quantum degen- 
eracy. The effect of phase space blocking significantly 
reduces the collision rates at quantum degeneracy. Time 
scales for evaporation are greatly prolonged for fermions, 
while for bosons the final state stimulation factor (l+n e ) 
further enhances the classical run-away evaporation. Of 
course this type of slowing is not unique to our scheme of 
fermionic evaporation, it is a generic feature of fermionic 
quantum statistics. As such, it is also present for other 
fermionic cooling schemes such as sympathetic cooling 



with bosons. One should realize that in this regime of 
phase space density the energy scale involved is already 
extremely low, and the crucial challenge is in keeping the 
inelastic collision and loss rates down. 

In numerical studies with Eq. (Q) for a variety for 
parameters, we have found that down to the limit of a 
trapped gas temperature of the order of the Fermi en- 
ergy up (typically ~ 0.5 fip), the evaporation proceeds 
as if the atoms were classical or bosonic, i.e., the fi- 
nal state blocking effect is unimportant. In Figs. |] 
and ^| we display comparisons of evaporative cooling for 
bosons, fermions, and classical atoms (with the same 
collision cross section a) contained in an isotropic har- 
monic trap with trap frequency hu>. The system con- 
tains N = 10 5 atoms initially, which are assumed to 
obey the Maxwell Boltzmann distributions with an ini- 
tial temperature T = 164.92 (Tuo) = 2Tp, i.e., twice 
the initial Fermi temperature. The evaporation is ef- 
fected by ramping down the trapping potential with a 
cut e cu t(t) = eoe~ 7cutt - The time scale of the evaporation 
is r = t/ro with 1/tq = (ma /ir 2 h 3 )(hui) 2 . The evapo- 
ration rate is 7 cu t = 7coii(0)/30, with 7 co ii(0) the initial 
collision rate. The time dependent collision rate is given 
by 



: 1 "» ( ' ' = mjS 12 ^Z+epeg+e/rie-, Cj, Cg, £f) 



-n e ,(t)n ej (t)[l+rin eR (t)}[l + V n er (t)}y (11) 

In Fig. [| the temperature Tit) is determined by fitting 
a Maxwell, Bose, or Fermi distribution function with a 
temperature T(t) and fugacity z to the total number of 
remaining atoms N(t) — ^ e p(e)n e [z(t), T(t)] and the 
average energy e = p(e)en e [z(f), T(f)]. In the trap, 
the Fermi temperature Tp and the critical temperature 
for Bose-Einstein condensation Tq are given by ksTp = 
Hlo(6N)^ 3 and k B T c = ?iw(iV/1.202) 1 /3. 



III. SYMPATHETIC COOLING OF HYPERFINE 
FERMIONIC MIXTURES 

In a recent experiment the JILA group has succeeded 
in implementing rf evaporation on a trapped fermionic 
cloud containing two spin states |l(|. The effect of the 
low energy p-wave suppression was clearly demonstrated 
in their measurements of the collision rates. Our gen- 
eral analysis of anisotropic collisions in this paper can 
be adapted to an analysis of this experiment, in par- 
ticular, by specializing to the case of s-wave collisions 
between the two trapped spin states, and p-wave colli- 
sions for each of the two spin states. We first discuss the 
appropriate collision integrals again assuming that trap 
confinement is not important during the collisions, i.e., 
when the ground state size of the trap is much larger 



than the effective range of the interatomic potential; this 
is true for almost all magnetic traps currently employed 
in laboratories. The scattering cross section 



\F B/F (q,cf)\ 2 dq' 

fdai ^ n dfa=*^dtf (12) 



V d v' /iab a 

is invariant in different reference frames, and this enables 
us to write the differential scattering cross section for p- 
wave scattering as 

= aq 4 cos 2 6a cos 2 6*, (13) 

dq 

where 9, 0' specify the directions of the relative momen- 
tum before and after the scattering. The constant a is 
obtained from the experimental data Jl(| in terms of an 
equivalent s-wave scattering cross-section <t p (Tq) — o s at 
a temperature To, i.e., 

0^(4%) =a s , (14) 

where q 4 is the average of q taken with respect to a 
Boltzmann distribution. This yields 

(7,1 , 

(15) 



15-7T (mfc^To) 2 

We can then substitute Eq. ( |l2| ) into (Q), and integrate 
over the phase space dfdpi again using the assumption 
of ergodicity 

f(f,p) = J den{e)8[e-p 2 /2m-U{i% (16) 

to obtain the ergodic Boltzmann equation for p-wave part 
of the scattering 

, . d , . mcr 8 1 f , . 

p{e A )-n{e A , t) = — J de, de 2 de 3 6(e 1 +e 2 -e 3 - e 4 ) 

P( e <) e ^3 (d ~ £2) 2 (£3 - 

n(ei)n(e a )[l + ??n(e 3 )][l + V»(e*)] - [1 + Vn[ei)][l + 7?n(e 2 )]n(e 3 )n(e 4 ) , (17) 

where e> (e<) is the largest (smallest) of the four ener- 
gies Cj. We see that the collisions become ineffective as 
soon as the energy scales become less than /c^Tq. To ob- 
tain a classical collision rate we substitute the Maxwell 
Boltzmann distribution 

n(e) = z- 1 e- e/kaT , (18) 

into equation |lj], set rj = 0, and replace the sums by 
integrals and integrate the first term on the r.h.s. over 
64. This yields the p-wave collision rate as 

, . TO(T, , . o IJT 
T„=#)-|N 2 ^. (19) 
?r 2 ft 45T 



with the scaled temperatures T = fcsT ' jhw and I p = 5 is 
the value of the fourfold integral 



dx\ dx2 dx% dx4<5(xi + x 2 — x 3 — X4)e e x2 



p(x<) 



3a;; 



(,; ; 



(Xl - X 2 ) 2 {X3 - XiY 



(20) 



In comparison the total collision rate for s-wave scatter- 
ing yields 



7, = ^v(t)^(M 2 =: 



(21) 



with I s = 0.5 is the value of the integral similar to ( pQ ) 
but only with p{x < ) in the second line. In the above all 
results were given for a spherically symmetric trap with 
trap frequency to. The case of anisotropic traps can be 
obtained by substitute u> with (ujxiUyUJz) 1 / 3 . 

On the other hand, the proposed p-wave collision (in- 
side the electric field) would not experience such a slow 
down. Inside a dc-electric field, the p-wave scattering 
has the same non-momentum (fc) dependence as for the 
s-wave. We will now detail the treatment for a single 
term Fp(q,<f) — in tlQY{ Q (q)Y 10 (q') in the homogeneous 
case. Using approaches similar to that lead to Eq. (17), 
we obtain 



, . d . . 
P( e 4j^n(e 4 ,i) 



96to 



h 3 (M 3 



J dei de 2 de 3 S(ei 



(ei- £2 ) 2 (e 3 -e4) 2 
(e>-6<) 3 



£2 — ^3 

dq- — 



-e 4 ) 
b 2 - 



T 



b 2 /q 2 ) 



2\3/2 



n(ei)n(e 2 )[l + r?n(e 3 )][l + r?n(e 4 )] - [1 



where b = + v / e 7 >)/(v / e> — >/£>)■ Although we have 
not been above to complete the inside integral analyti- 
cally, we can immediate see that the energy dependence 
for the collision rate is such that a run away evaporation 
is possible in the classical limit, i.e. 



1 + b 2 - q 2 
r)n(ei)][l + J7n(e 2 )]n(e3)n(e4) 



. (22) 



ip 



(M =■ 



y ' ir 2 h 3 



(23) 



We have simulated the evaporation process for the re- 
cent experiment at JILA jl6| . We assume that two hyper- 
fine states were in the same magnetic trap, and include 
two types of collisions: the s-wave collision between the 
two states and the p-wave collision within each hyperfine 
state as discussed above. Inelastic processes are ignored. 
We start with 10 7 particles in each state at an initial 
temperature which is 135 times the initial Fermi temper- 
ature. Then we evaporatively cool the mixture by cut- 
ting off the hot atoms (from both states) with a cut rate 
which is 1/20 times the initial collision rate. The atoms 
rethermalize via p-wave and s-wave scattering. The p- 
wave collision rate decreases as the temperature of the 
gas decreases whereas the s-wave collision rate increases 
due to run-away evaporative cooling down to a temper- 
ature when Fermi-statistics become important and the 



final state blocking term 1 — n(e) causes the collision rate 
to decrease. At the end of the cooling process the temper- 
ature is more than 9 times smaller than the final Fermi 
temperature and there are less than 1.3 x 10 6 particles 
left in each state. This corresponds to an increase of 
phase-space density of more than 10 12 . The final average 
peak density in the trap is then n = 4.0 x 10 12 /cm 3 

In Fig. H we show our results for the numerically com- 
puted collision rates. We note that the p-wave colli- 
sion rate essentially dies out at the temperature mea- 
sured by the JILA experiment jl6) . The sum of the two 
scattering rates initially decreases due to the suppressed 
p-wave collisions, the total rates eventually increases as 
the runaway evaporation sets in for the s-wave scatter- 
ing. We started the evaporation at an initial temper- 
ature of Ti = 5.23 x 10 5 and have used the parame- 
ters of a s = a_p = 10~ n cm 2 at T Q = T;/4 from the 
JILA results [n6|. Assuming an average trap frequency 
of uj as (2tt) 80 (Hz) this corresponds to an initial temper- 
ature of Ti = 200 (pK) . The simulation was started with 
500 energy bins and each bin covers an energy interval of 
800 Tuo. After the cut-energy is ramped down to 40 en- 
ergy bins, the energy interval is decreased by a factor 10, 
a redistribution of all particle numbers to 500 finer energy 
bins is required. This procedure is repeated twice before 
one finally gets down to the required low temperature. 
The accuracy and convergence were checked by simula- 
tions with larger number of energy bins. In Fig. || we 
show the temperature of the evaporated mixture as well 
as the corresponding Fermi temperature [computed in 
terms of N(t) and T(t)]. We note that with appropriate 
initial conditions the gas will cross the Fermi degeneracy, 
at a time idcg = 3.3 x 10~ 3 u; 2 for the parameters used in 
the simulation , e.g. 827 (s) for u>/2ir = 80 Hz. We have 
not included inelastic loss processes such as three body 
collisions, however as long as the inelastic loss rate G is 
such that Gn 2 < , the cooling process will success- 
fully reach the degenerate regime. 



IV. CONCLUSION 



In conclusion, we have proposed the possibility of evap- 
orative cooling of trapped fermionic atoms in an external 
electric field. We have derived effective quantum kinetic 
equations for fermionic atoms with anisotropic dipole- 
dipole interactions, and presented numerical simulations 
for fermionic evaporative cooling. We have also discussed 
sympathetic cooling of fermionic mixtures, related to a 
recent experiment [ p.6| , when two hyperfine states of 40 K 
are magnetically trapped. Cooling is principally due to 
the predominant s-wave collisions between the two states. 
We have observed intrastate p-wave suppression consis- 
tent with the experimental results. In the region where 
degeneracy becomes important we have contrasted the 
qualitative differences in the evaporative cooling dynam- 
ics of bosons and fcrmions. 
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FIG. 1. The low energy collision cross section for spin 
triplet 6 Li (fcrmion). The dashed horizontal line indicates 
the field free value of triplet 7 Li [boson, a BC = —27.6 (a.u.)]. 
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FIG. 2. The same as in Fig. 1 but for 40 K (fermion). 
The dotted line is the field-free reference value for 39 K [boson 
a sc ~ 80 (a.u.)] and the dashed line the field-free reference 
value for 41 K [boson a BC ~ 300 (a.u.)]. 
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FIG. 3. The same as in Fig. 1 but for Rb (fermion) which 
is similar to the results of 82 Rb and 84 Rb. The dotted line 
is the field-free reference value for 85 Rb [boson, a BC ~ 369 
(a.u.)] and the dashed line the field-free reference value for 
87 Rb [boson, a sc ~ 106 (a.u.)]. 
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FIG. 4. Comparisons of the evaporative cooling for bosons, 
fermions, and classical atoms. The remaining number of 
atoms in the trap. The steps in the number of fermions for 
larger times is due to cutting discrete energies for a distribu- 
tion function which does not rethermalize fast enough. 




FIG. 5. The same as in Fig. 4, but for the numerically 
computed collision rate. 




FIG. 6. The same as in Fig. 4, but for the fitted tempera- 
tures of the remaining atomic distributions. 




FIG. 7. The same as in Fig. 4, but for the average energy. 




FIG. 8. The numerically computed collision rates j p and 
7 S for the discussed two state evaporation simulation. 




FIG. 9. The temperature and the Fermi temperature for 
the evaporated two state Fermi mixture. 



